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ABSTRACT 


Some  elementary  problems  peculiar  to  a 
ground-to-ground  homing  missile  system  are 
discussed.  The  work  is  specifically  oriented 
toward  the  application  of  a  television  guidance 
sensor  of  the  contrast  contour  type  to  an  anti¬ 
tank  weapon.  The  design  of  a  simple  laboratory 
TV  servo-tracker  is  presented.  Appendices 
A  and  B  discuss  test  vehicles  for  subject 
applications  and  present  pertinent  XQ-2  drone 
parameters. 
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SOME  OPTICAL  CONTRAST  SEEKER 
SYSTEM  CONSIDERATIONS 


I.  INTRODUCTION 

The  employment  of  passive  seekers  of  the  television  (raster-scan) 
type  as  the  error  sensing  element  in  a  collision-seeking  (homing)  missile 
presents  many  peculiar  problems.  A  thorough  investigation  of  the  system 
application  of  such  a  seeker  dictates  the  need  for  specific  information 
in  the  following  areas: 

1.  Acquisition  and  identification  of  targets. 

2.  Optimum  selection  of  optics  and  sensor  and  evaluatiox.  of  the 
limitations  imposed  by  the  operating  environment. 

3.  Signal  processing  techniques  for  acquisition,  track,  stability 
and  communication  (if  needed). 

4.  Practical  seeker  design  specifications  to  meet  dynamic  require¬ 
ments  and  impose  minimum  size  and  weight  requirements  on  the  missile. 

5.  Missile  navigation  and  flight  dynamics  which  relate  to  seeker 
requirements . 

This  report  covers  the  work  done  to  date  on  areas  4  and  5.  Addi¬ 
tional  related  work  has  been  reported  in  Report  No.  RE-TR-63-3. 

For  purposes  of  limiting  the  number  of  system  parameters,  a  fixed- 
range  ground-to-ground,  line -of-sight  (LOS)  tank  engagement  was 
assumed.  Within  this  framework,  simplified  mathematical  missile 
models  were  used  in  a  proportional  navigation  homing  environment  for 
the  establishment  of  gross  system  requirements,  limitations  and 
capabilities.  The  design  data  thus  generated  were  used  to  outline  the 
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requirements  of  a  two-axis  servo  tracker,  and  a  laboratory  model  of 
the  tracker  has  been  constructed. 

II.  SOME  ELEMENTARY  GEOMETRIC  CONSIDERATIONS 


In  the  design  of  a  homing  system  based  on  an  image  contour  tracker 
of  the  raster-scan  type,  the  angle  subtended  by  the  target  within  the 
field-of-view  becomes  interesting  for  several  reasons: 

1.  At  maximum  range,  the  field-of-view  of  the  optics  must  be 
sufficiently  small  to  yield  a  target  image  large  enough  to  exceed  the 
resolution  limitations  of  the  TV  seeker. 


2.  As  the  missile  nears  the  target,  the  image  will  subtend  a 
proportionally  larger  portion  of  the  field-of-view,  until  eventually  the 
image -to -field  area  ratio  becomes  unity,  and  tracking  information  is 
lost.  The  missile  must  then  coast  unguided  to  the  target.  Minimizing 
this  coast  range  requires  a  large  field-of-view  conflicting  with 
resolution  requirements. 


This  conflict  may  be  minimized  by  using  programed  variable -field 
optics  at  a  penalty  in  cost,  weight,  complexity,  and  size.  To  establish 
allowable  terminal  engagement  guidance  system  performance ,  the 
trade-offs  involved  in  using  fixed  field  optics  were  investigated.  For 
this  purpose,  the  magnitude  of  the  angle  subtended  by  a  given  target  as 
a  function  of  range  is  required  and  was  determined  as  follows: 


From  Figure  1: 


et  =  2\ 


=  tan 


-l  h/2 
Rt 


Since  X.J-  is  small  for  all  ranges  of  interest,  the  usual  small-angle 
approximations  apply.  For  example,  if  h/2  =  5  ft.  ,  Xt  will  equal 
100  milliradians  or  6”  at  R  =  50  feet. 


Thus: 


and 
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It  can  be  seen  that  is  dependent  only  on  target  size  and  range.  Since 
the  above  function  is  a  reciprocal  function,  it  is  not  conveniently  plotted 
on  linear  graph  paper.  However,  one  can  make  the  transformation: 

l°g  6t  ~  log  l1  -  log  R  . 

Ignoring  the  logarithm,  the  above  equation  is  obviously  that  of  a 
straight-line  (slope-intercept  form,  Y  =  b  -  mx).  Thus,  if  the  .function 
0t  =  h/R  is  plotted  on  log-log  paper,  straight-line  plots  are  obtained. 
The  advantages  of  this  form  for  the  given  function  are: 

1.  Several  decades  of  the  variables  of  interest  are  easily  plotted 
with  good  accuracy. 

2.  Only  three  values  of  each  function  of  interest  need  be  calculated 
to  verify  the  straight  line.  Plots  of  versus  R  are  shown  in  Figure  2 
for  three  different  target  sizes.  Note  that  four  decades  are  used  for 
the  ordinate.  The  family  of  curves  to  the  right  should  be  referred  to 
the  right  hand  ordinate,  the  one  to  the  left,  to  the  left  hand  ordinate. 

For  the  range  of  fields  of  view  considered  appropriate  (0.5°  to  2°), 
Figure  2  shows  that  the  coast  range  could  be  as  great  as  1, 000  feet  if 
fixed  optics  are  used  with  the  proposed  TV  tracker--  a  contrast  contour 
device  (CCT  tracker)  being  built  for  the  Army  by  the  Southern  Research 
Institute.  Terminal  performance  of  missile  models  were  evaluated 
under  such  conditions,  and  some  preliminary  results  are  reported  in  a 
later  section. 

In  addition  to  6^,  it  might  be  of  interest  to  know  the  rate  of  change 
of  the  angle  subtended  by  the  target,  or  0j.  This  can  be  easily 
determined: 

Since,  6t  = 

•  d  /h  \  hR  hv  Ov 
9t  ~  dt  \R/  ~  rT  ""  R2  °r  R  ' 

Unlike  0).,  is  dependent  upon  missile  velocity.  Figures  3  and  4  show 
plots  of  for  three  target  sizes  and  two  missile  velocities  (750  ft/ sec 
and  1,000  ft/sec,  respectively).  For  the  reasons  previously  indicated, 
these  plots  are  also  made  on  log-log  paper. 

Missile  velocity  was  assumed  constant  throughout  each  flight  for 
all  analyses  reported. 
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Another  parameter  of  interest  is  the  angular  turning  rate  required 
of  a  missile  in  order  to  track  a  moving  target.  Since  the  maximum  rate 
corresponds  to  a  target  moving  perpendicularly  to  the  missile  flight 
path,  approximate  values  were  obtained  thus:  At  any  one  instant,  the 
target  was  considered  to  be  a  point  moving  along  the  circumference  of 
a  circle  of  radius  R  (Fig.  5).  The  target  velocity  was  then  V^,  the 
tangential  velocity  of  the  point,  and 

Vt 

angular  velocity  of  missile  ~  =  -r~  , 


angular  acceleration  of  missile 


R1 2 3 4  R2 


co  versus  R  is  plotted  in  Figure  6  and  a  versus  R  in  Figure  7.  In  these 
plots,  target  velocities  of  44  ft/sec  and  88  ft/sec  are  considered.  For 
the  30  mph  (44  ft/sec)  target,  moving  perpendicularly  to  the  field  of 
view,  the  required  tracking  rates  are  quite  modest  to  within  a  very 
short  rang'e  from  the  target.  Actual  angular  rates  will  include  effects 
of  missile  dynamics,  which  will  also  affect  seeker  angular  rate 
performance  requirements. 


III.  MISSILE  STUDIES 


A.  Introduction 


The  objective  of  the  present  Optical  Contrast  Seeker  (OCS) 
program  is  to  demonstrate  the  feasibility  of  a  homing  missile  system 
employing  such  a  seeker,  rather  than  to  develop  a  tactical  system 
design.  Many  simplifying  assumptions  are  therefore  possible  and 
indeed  desirable.  Preliminary  analytical  and  analog  simulation  studies 
were  made  within  the  framework  of  the  following  limitations: 

1.  Missile  dynamics  are  represented  by  a  second  order  transfer 
function. 

2.  Planar  analysis  only  is  performed. 

3.  Target  engagement  is  line -of  -  sight  from  a  stationary  launcher. 

4.  Target  is  stationary. 
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Homing  missiles  make  use  of  the  observed  motion  of  the  missile- 
target  LOS,  and  several  well-known  navigation  laws  are  available  for 
selection.  Among  these  are  pure  pursuit,  constant  bearing,  and  pro¬ 
portional  navigation.  The  first  two  generally  require  extreme  missile 
maneuvers;  for  vehicles  of  limited  performance  the  navigation  scheme 
is  the  most  practical  to  implement  and  is  thus  chosen  here  (Ref.  1).  In 
proportional  navigation,  the  missile  detects  any  rotation  of  the  line  of 
sight  and  assumes  a  turning  rate  proportional  to  the  line  of  sight  rotation, 
in  the  direction  which  will  reduce  the  LOS  rate  and  thus  tend  to  return 
to  a  constant-bearing  collison  course.  In  the  case  of  a  nonaccelerating 
target,  once  the  missile  has  corrected  initial  errors,  the  line-of-sight 
ceases  to  rotate  and  the  constant  bearing  course  is  achieved. 

Ideally,  proportional  navigation  can  be  defined  thus:  =  Kcr , 

i.  e.  ,  the  missile  turning  rate  is  equal  to  some  constant  (usually  between 
2  and  6)  times  the  LOS  angular  rate.  A  high  value  of  K  relates  to  loop 
gain  and  generally  yields  an  oscillatory  trajectory.  On  the  other  hand 
for  K  =  2,  it  can  be  shown  mathematically  that  the  required  missile 
acceleration  is  no  longer  bounded  and  approaches  infinity  as  the  range 
approaches  zero  (Ref.  1). 


If  a  second  order  missile  and  first  order  lag  for  the  seeker  are 
assumed,  the  guidance  equation  would  become: 


'Sz  |  2^ 


1  +  T  ,S 


This  equation  can  be  derived  from  a  differential  equation  analysis: 
Assumed  approximate  missile  behavior: 

a  +  2£«n  a  +wn2a  = 

where  relates  to  missile  control  effectiveness.  Replacing  d/dt  by 
S  and  dividing  through  by  wnz  : 


+  —  s  +  i 


wn 

Assuming  the  seeker  equation  with  1st  order  lag: 
Tj  6  +  6  =  Kir 
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6(1  +  TlS)  =  K<rS  ,  or  6-  ■  K— S--  . 

1  +  Ti  b 

Inserting  6  into  Equation  A  and  transposing  terms  in' order  to  obtain 
the  Oi  /<r  transfer  function: 


Introducing  the  approximate  lift  equation: 

.  O' 

r2 


=  ~  w^ere  T2  assumed  constant,  or  or  - 


Inserting  above  value  of  a  into  Equation  B  the  desired  y^/cr  transfer 
function  is  obtained: 


_M 

tr 


KKj 

t2 


SL+2is  +  J  I1  +  rxSl 


GJr 


The  above  expression  has  the  same  form  as  the  original  equation. 
By  transposing  u  to  the  right  side  and  combining  it  with  the  free  S,  crS 
or  <r  is  obtained. 


if 

The  block  diagram  of  this  simple  homing  loop  is  shown  in  Figure  8, 
and  the  analog  simulation  program  for  one  plane  (yaw)  in  Figure  9. 
Several  elementary  observations  can  be  made  from  Figure  8: 


1.  The  system  is  of  the  4th  order,  and  thus  the  loop  will  inevitably 
become  unstable  if  the  gain  becomes  high  enough.  This  is  due  to  the 
fact  that  the  eventual  contribution  of  a  4th  order  system  will  be  360° 
of  phase  lag.  However,  provided  l/ij  is  greater  than  the  natural 
frequency  of  the  missile,  the  total  phase  lag  in  some  region  below  u>n 
is  less  than  180°,  and  the  system  will  operate  stably  in  this  region. 


Z.  Assuming  that  small  angle  approximations  are  valid,  <r  is 
given  by  y/R.  This  means  that  as  R  approaches  0,  y/R  approaches 
00  and  the  system  will  of  necessity  become  unstable.  Thus,  for  a 
given  missile  performance,  the  choice  of  navigation  constant  K 
becomes  a  matter  of  trade-offs  between  launch  and  terminal  performance. 
The  initial  gain  must  be  high  enough  so  that  dispersion  is  minimized, 
yet  low  enough  so  that  the  system  does  not  become  unstable  prematurely. 

3.  The  navigation  constant  K  tends  to  lose  explicit  meaning  in  this 
type  of  simulation.  All  the  constants  K,  K1(  t2,  are  lumped  together 
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as  a  factor  relating  to  the  loop  gain.  Only  when  specific  performance 
numbers  are  assigned  to  the  missile  (K^,  72).  does  K  become  meaningful. 
For  convenience  in  mechanizing  the  simulation,  K  was  assumed  to  be 
3,  and  K1(/t2  (X.  in  Figure  9)  was  varied. 

C.  Investigation  of  Fixed  Versus  Gimballed  Seeker 

An  early  simulation  study  was  required  in  order  that  an  expedi¬ 
tious  decision  could  be  made  as  to  whether  a  constant  field  of  view 
seeker  could  be  body  fixed.  The  dominant  requirement  of  a  body-fixed 
tracker  is  that  the  missile  center-line  (or  other  longitudinal  reference) 
form  an  angle  with  the  LOS  not  greater  than  the  field  of  view  throughout 
the  flight;  otherwise,  the  target;  will  leave  the  field  of  view  and  tracking 
will  terminate.  At  the  time  this  information  was  desired,  it  was  believed 
that  a  field  of  view  of  ±  2°  might  be  feasible  and  the  studies  were  per¬ 
formed  on  this  basis.  A  much  smaller  field  of  view  is  now  required  to 
allow  the  CCT  tracker  to  resolve  tank-sized  targets  at  maximum  range. 
However,  it  was  determined  that  even  a  ±  2°  field  of  view  would  not 
allow  employment  of  a  body-fixed  seeker  unless  some  stringent  missile 
launch  dispersion  restrictions  were  met.  For  the  smaller  field  of  view, 
a  servo-directed  seeker  sensor  is  a  definite  requirement. 

For  the  purpose  of  investigating  launch  dynamics,  a  series  of 
simulated  trajectories  were  obtained  for  the  yaw  plane.  Figure  10, 
showing  pertinent  geometry,  indicates  that  the  angle  a  +  y  is  of  the 
greatest  interest.  (3((3  =  o')  can  be  neglected  since  it  is  generally  small 
and  its  peak  value  occurs  after  O'  and  y  have  been  reduced  considerably. 

The  fixed  parameters  used  in  the  runs  were; 

R0  =  initial  range  =  6,  000  feet 
YQ  =  initial  cross  velocity  =  50  ft/sec 
Tj  =  seeker  time  constant  =.  0.  03  sec 
K  =  navigation  constant  =  3 
Stationary  target 

The  initial  cross-velocity  of  50  ft/sec  was  selected  to  take  care  of 
all  dispersions  that  might  be  encountered.  This  value  is  less  than 
5  per  cent  of  the  missile  velocities  selected  and  is  considered  to  be  a 
realistic  estimate.  All  other  parameters  were  varied  and  are  tabultated 
in  Table  I.  It  was  found  from  the  simulator  runs  that  the  Y0  used 
resulted  in  an  initial  y  that  did  not  change  noticeably  for  a  few  hundred 
feet.  The  peak  transient  value  of  a  occurred  generally  at  T  <0.3  seconds; 
thus  the  maximum  angle  between  missile  centerline  and  target  can  be 
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approximated  by  (a  +  yQ  +  |3)  ~  (a  +  yo).  Table  I  shows  that  without 
exception,  all  angles  exceed  the  2°  half-angle  of  the  assumed  4°  body- 
fixed  seeker. 

No  analysis  was  made  of  terminal  guidance  in  this  study.  Guidance 
was  terminated  and  the  controls  locked  when  the  missile  was  50  feet 
from  target.  This  distance  was  arbitrarily  selected  and  for  this  reason 
many  of  the  runs  resulted  in  quite  unstable  terminal  trajectories.  Thus 
the  miss -distance  column  of  Table  I  is  not  too  meaningful,  and  for 
those  runs  where  the  oscillations  were  very  large  (i.e.  ,  loop  became 
quite  unstable)  prior  to  the  50-foot  range,  no  values  were  entered.  The 
basic  reason  for  omitting  terminal  guidance  in  this  study  is  based  on 
the  fact  that  for  small  fields  of  view  tracking  information  will  be  lost 
in  the  terminal  phase  and  the  missile  will  coast  without  guidance  as 
previously  discussed. 

As  a  result  of  these  simple  investigations,  some  preliminary 
conclusions  may  be  hazarded: 

1.  Unless  a  fixed  seeker  with  a  half-angle  of  at  least  5°  to  7°  is 
feasible,  further  consideration  should  be  given  only  to  gimballed 
seeker  sensors.  Subsequent  work  is  based  on  this  conclusion. 

2.  Launch  dynamics  place  performance  requirements  that  are  not 
unreasonable  on  the  servo  tracker.  A  typical  peak  value  of  2°  for  an 
underdamped  (£,  =  0.  2)  missile  may  be  reached  in  a  maximum  of  0.  1  sec. 
An  approximate  ot  of  20°/sec  then  results,  and  we  must  track  at  least 

at  that  rate  to  stay  on  target.  At  this  time  it  seems  that  a  gross  mini¬ 
mum  &  value  of  30°/sec  would  be  adequate  as  a  specification  provided 
further  guidance  studies  do  not  dictate  otherwise.  A  more  accurate 
requirement  can  be  determined  only  for  well  defined  missile  performance 
parameters.  Higher  rates  up  to  60°/sec  should  be  achievable  without 
undue  difficulty  if  necessary. 

D.  Stability  Considerations 

A  detailed  stability  analysis  of  any  feedback  control  system 
requires  specific  dynamic  performance  of  each  element  of  the  control 
loop,  including  nonlinearities.  The  determination  of  the  degree  of 
stability  of  the  homing  loop  shown  in  Figure  8  as  a  function  of  range 
can  be  determined  easily  but  should  not  be  construed  as  indicating  the 
performance  of  an  actual  missile  system.  Such  an  analysis  at  this 
stage  is  nevertheless  valuable  for  several  reasons: 
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Table  I.  LAUNCH  DYNAMICS  DATA 


I 


! 


1 

I 

I 


1.  Stability  analysis  of  an  idealized  linear  system  generally  rep¬ 
resents  what  could  be  called  a  "best"  performance,  and  a  similar 
physical  system  could  be  expected  to  have  poorer  response.  Thus  a 
broad  upper  bound  of  performance  can  be  established. 

2.  The  postulate  that  tracking  information  will  probably  be  lost 
with  the  OCS  tracker  loop  before  loop  stability  becomes  a  problem 
can  be  verified  by  correlating  the  critical  loop  gain  in  terms  of  range 
with  the  information  shown  in  Figure  2. 

The  overall  open  loop  transfer  characteristics  (HG(S))  of  the  system 
shown  in  Figure  8  are  given  by: 

KKtVafr2  S 

HG(S)  =  (S2  +  2 1,  G)n  S  +  wn2  )(l  +  t i  S)  S2 
The  value  of  loop  gain  A  is  given  by: 

A 

t2R 

If  the  value  of  A  at  which  the  system  goes  absolutely  unstable  is 
determined,  the  range  at  which  this  occurs  is  easily  computed  by  means 
of  the  above  equation  since  all  other  terms  are  known  (assumed). 

The  Routh  Stability  Criterion  (Ref.  2)  is  probably  the  easiest  means 
of  establishing  the  value  of  gain  for  critical  stability.  However,  if  the 
root  locus  plots  of  the  various  assumed  systems  are  obtained,  one  can 
not  only  obtain  the  critical  gain,  but  also  values  of  gain  (and  consequently 
R)  in  the  stable  region.  For  example,  it  might  be  of  interest  to  deter¬ 
mine  the  range  at  which  a  damping  ratio  of  0.  1  occurs.  For  this  reason, 
the  root  locus  method  was  used  in  this  study.  Values  of  u)n  and  t,  were 
chosen,  representing  four  different  missile  models;  a  seeker  lag  time 
constant,  navigation  ratio,  and  missile  performance  parameters  were 
assumed  for  a  1,  000-ft/sec  missile.  The  root  loci  of  these  systems 
are  shown  in  Figures  11  through  14,  and  the  values  of  R  for  various 
degrees  of  stability  are  shown  in  Table  II.  Table  III  summarizes  a 
comparison  of  the  critical  range  Rc  with  the  value  of  R  at  which  a 
7/2  X  7/4 -foot  target  image  fills  various  angular  fields  of  view.  This 
elementary  analysis  tends  to  support  the  postulate  that  the  sensor  may 
very  well  be  the  limiting  element  in  terminal  guidance. 


10 


i 


Table  II 

STABILITY  VERSUS  RANGE 


o 

i-H 

II 

G 

3 

con  =  20 

Ri(ft) 

%imit 

R^  =  0.1 

R^  =  0.  2 

£  =  0.  1 

1, 350 

2,  780 

4  =  0.2 

834 

1,610 

2,  780 

£  =  0.1 

52.  4 

115 

£  =  0.2 

36.  4 

70 

112 

Table  III 

CRITICAL  RANGE  VERSUS  FIELD  OF  VIEW  FOR  TYPICAL  MISSILES 


0t 

degree 

Coast 

range, 

ft 

Critical  range 

wn  = 

20 

a>n  - 

10 

tr* 

II 

o 

i— 

£  =  0.  2 

\T« 

(1 

O 

H- ' 

£=0.2 

0.  25 

1,720 

52.  4 

36. 4 

1, 350 

834 

0.  5 

860 

52.  4 

36. 4 

1,350 

834 

1.0 

430 

52.  4 

36. 4 

1,350 

834 

2.  0 

215 

52.  4 

36. 4 

1,350 

834 

E.  Terminal  Phase  Studies 

Since  the  likelihood  that  the  missile  will  have  to  coast  into  the 
target  for  some  significant  range  had  been  established,  the  effects  on 
accuracy  of  this  unguided  trajectory  was  interesting.  For  simplicity, 
the  determination  of  miss -distance  was  broken  down  into  two  groups: 
stationary  targets  and  moving  targets. 

1.  Stationary  Targets.  The  definition  of  proportional  naviga¬ 
tion  implies  that  detection  of  the  line -of  -  sight  rotation  is  performed,  and 
this  signal  is  used  to  correct  the  missile  flight  path,  the  tendency  being 
to  approximate  the  constant  bearing  course.  Therefore,  in  the  case  of 
a  stationary  target,  once  initial  heading  errors  are  corrected,  the  line- 
of-sight  ceases  to  rotate  and  a  constant  bearing  course  is  achieved.  If 
cr  is  zero,  then  obviously  y  is  also  zero,  and  the  missile  will  remain 
on  the  correct  collision  course  throughout  the  flight  (assuming  no 
external  disturbances),  since  no  guidance  information  is  needed  once  a 
steady  state  cr  value  is  achieved.  Thus  the  displacement  of  the  initial 
LOS  decreases  as  a  function  of  time-to-go,  becoming  zero  at  T  =  0. 
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The  obvious  limitation  on  accuracy  during  the  coasting  period  on  a 
stationary  target  is  the  fact  that  no  error  is  sensed  if  the  missile  devi¬ 
ates  from  its  correct  trajectory  due  to  disturbances;  these  can  be 
broken  down  into  two  distinct  types:  random  and  drifts.  Due  to  the 
principle  of  superposition  applicable  to  linear  systems,  the  effect  of 
each  type  of  disturbance  may  be  evaluated  individually  and  the  total 
effect  is  obtained  if  desired  by  the  linear  combination  of  the  individual 
ones. 

An  evaluation  of  the  effects  of  random  disturbances  was  performed 
by  superimposing  upon  the  error  signal  e  a  random  noise  signal  with  a 
Gaussian  distribution  about  0.  Three  amplitude  values  were  selected 
to  correspond  to  various  reasonable  signal  to  noise  ratios.  Since  for  a 
stationary  target  cr  becomes  a  constant,  the  summing  of  a  constant  rms 
amplitude  noise  signal  yields  a  fixed  signal  to  noise  ratio. 

At  the  moment  guidance  is  terminated,  the  missile  heading  will  be 
given  by  cr  ±  e  where  cr  is  the  "correct"  heading  and  6  is  the  error  caused 
by  the  instantaneous  noise  amplitude  at  the  time  of  guidance  termination. 
The  effect  of  the  heading  deviation  may  be  assessed  with  the  help  of 
Figure  15.  The  maximum  amplitude  that  e  may  have  at  the  moment 
guidance  is  terminated  and  a  target  hit  is  still  obtainable: 

h/2 

where 

Y  =  half  of  target  height 
Rc  =  coast  range 

Since  the  peak  amplitude  and  Gaussian  distribution  of  the  injected  noise 
are  known,  the  cumulative  probability  distribution  of  the  noise  signal 
may  be  plotted  as  shown  in  Figures  16  through  18.  These  plots  are  for 
the  three  selected  values  of  S/N  ratios.  The  percentage  of  time  that  € 
is  less  than  the  maximum  permissible  value  can  then  be  read  off  these 
curves  and  an  approximation  to  the  hit  probability  can  be  obtained. 

The  primary  limitation  of  this  estimate  is  that  a  steady  state  con¬ 
dition  for  the  controlled  system  is  assumed.  Although,  this  is  true 
with  respect  to  the  basic  error  signal  (cr  =  constant  at  this  stage),  the 
system  will  be  in  a  transient  state  due  to  the  noise,  and  the  actual  output 
(heading)  error  could  exceed  the  input  considerably  due  tb~the\  damping 
ratio  decrease  as  a  function  of  lessening  range.  However,  once  again 
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an  upper  bound  for  performance  is  obtained.  The  predicted  hit  proba¬ 
bilities  for  three  different  coast  ranges  and  signal-to-noise  ratios  are 
tabulated  in  Table  IV,  along  with  experimental  results  obtained  as 
follows:  In  order  to  obtain  an  idea  of  the  accuracy  of  this  prediction, 
an  analog  simulation  experiment  was  performed.  Approximately 
1, 000  trajectories  were  made  for  each  of  the  conditions  in  Table  V  and 
the  total  runs  and  hits  counted.  The  ratio  of  the  two  yields  the  hit  proba¬ 
bility.  A  hit  was  defined  as  any  trajectory  within  ±  3. 5  feet  of  target 
at  R  =  0.  A  sample  recording  of  part  of  one  of  the  runs  is  shown  in 
Figure  19.  With  the  exception  of  one  obvious  discrepancy,  Table  IV 
shows  that  fair  agreement  exists  between  predicted  and  experimental 
value  s . 


Table  IV 

COMPARISON  OF  EXPERIMENTAL  Ph 
VERSUS  ANALYTICAL  PREDICTION 


<r  =  28  mr 

R-coast<  ft 

Predicted  P^ 

Measured  P 

1 

S/N  =  7 

S/N  =3.5 

S/N  =  2 

S/N  =  7 

S/N  =3.5 

S/N  =  2 

800 

0.  95 

0.  6 

0.  3 

0.  88 

0.  7 

0.6 

400 

1.0 

0.9 

0.  7 

0.  87 

0.  76 

0.  65 

ZOO 

1.  0 

1.0 

0.  95 

0.9 

0.87 

0.  8 

Table  V 

MAXIMUM  ALLOWED  e  VERSUS  COAST  RANGE 


^■coast’ 

Maximum  allowable  e 

800 

4.  68  mr 

400 

9-  37 

200 

18.  74 

Unilateral  drifts  appear  to  present  the  most  serious  problem  during 
coasting  flight.  A  specification  of  maximum  allowable  cumulative  drift 
error  is  easily  obtained.  The  external  drift  must  impart  upon  the 
missile  a  lateral  acceleration  no  greater  than  that  which  will  result  in 
a  lateral  displacement  of  3.  75  feet  at  R  =  0.  This  value  is  given  by  the 
following  equation: 
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for  a  standard  7.  5  X  7.  5 -foot  target. 


2yV2  7.5  V2 

ad  =  R7"  =  Rc* 

The  limit  values  of  ad  for  a  1,  000  ft/sec -missile  velocity  and  three 
coast  ranges  are  shown  in  Table  VI. 


Table  VI 

MAXIMUM  ALLOWABLE  ad  VERSUS  COAST  RANGE 


^■coast>  ft 

a^,  ft/sec2 

800 

11.7 

400 

46.  8 

200 

182.  0 

2.  Moving  Targets.  The  problem  of  hitting  moving  or  accele¬ 
rating  targets  presents  the  most  serious  limitation  on  a  weapon  that  loses 
guidance  information  in  the  closing  phases  of  engagement. 

In  addition  to  noise  and  drift  errors,  the  error  caused  by  the  change 
in  target  position  during  the  coast  period  must  be  added.  The  magnitude 
of  this  error  is  a  function  of  many  parameters,  and  since  the  generalized 
proportional  navigation  trajectory  equations  cannot  be  solved  in  closed 
form,  a  general  expression  for  miss  becomes  difficult  to  obtain.  Some 
simulation  trajectories  for  the  typical  missile  transfer  functions  pre¬ 
viously  used  have  indicated  that  miss  distances  of  1  0  feet  are  obtained 
against  accelerating  (4.5  ft/sec2)  targets  for  an  800-foot  coast  range. 
Lesser  target  accelerations  and  shorter  coast  ranges  will  diminish  this 
error  proportionally. 

Constant  velocity  targets  present  the  least  problem,  since  the 
dynamics  involved  are  very  similar  to  an  initial  pointing  error.  Most 
of  the  correction  is  performed  early  in  flight  (for  a  reasonable  loop 
gain)  and  the  miss-distance  due  to  lack  of  terminal  guidance  is  minimized. 

In  summarizing,  it  may  be  stated  that  the  overwhelming  control 
problem  affecting  the  success  of  employing  a  TV  contour  tracker  is 
guidance  loss.  Against  stationary  targets  this  problem  could  conceiv¬ 
ably  be  circumvented  by  a  refined  missile  and  guidance  system  design 
which  would  minimize  dirft  and  noise  disturbances.  Successful  engage¬ 
ment  of  accelerating  targets  however  would  seem  to  require  that  the 
coast  range  be  minimized,  ideally  to  negligible  distances.  An  exact 
quantitative  determination  of  this  probability  against  accelerating 
targets  as  a  function  of  miss -distance  would  require  a  much  more 
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elaborate  analog  mechanization  than  the  present  simulation  and  will 
form  part  of  future  work. 

It  should  be  pointed  out,  lest  previous  statements  leave  an  errone¬ 
ous  impression,  that  loop  instability  need  not  be  a  limitation  on  terminal 
guidance.  Modifications  of  the  proportional  navigation  scheme  are 
possible  by  the  judicious  tailoring  of  loop  gain  to  allow  accurate  guid¬ 
ance  essentially  to  the  impact  point.  This  has  been  done  in  the  past, 
particularly  in  the  air-to-air  missile  field  where  the  dynamic  problem 
is  much  more  severe.  The  usual  approach  has  been  to  design 
the  missile  and  control  loop  around  carefully  defined  engagement  require¬ 
ments.  This  of  course  is  necessary  within  the  limits  of  physical 
realizability.  Additionally,  for  a  missile  with  limited  lateral  accelera¬ 
tion,  instability  may  result  in  limit-cycle  oscillation  rather  than  in  a 
diverging  trajectory;  if  the  natural  frequency  is  sufficiently  high  the 
amplitude  of  these  oscillations  can  be  made  quite  small.  Thus  instability 
need  not  necessarily  result  in  misses. 

IV.  SERVO  TRACKER  DESIGN 

Since  the  high  probability  has  been  established  that  a  seeker  system 
utilizing  the  CCT  type  of  tracker  will  require  a  servo  gimbal  mount, 
an  experimental  device  was  assembled.  An  existing  gimbal  system 
was  adapted  for  use  as  a  laboratory  model  of  a  TV  tracker.  The  CCT 
equipment  was  not  available  for  this  purpose;  however,  a  Walleye 
type  tracker  was  substituted.  This  is  a  point  contrast  tracker  which 
derives  positional  information  (with  respect  to  an  arbitrary  reference) 
by  means  of  slewing  vertical  and  horizontal  gates,  the  intersection  of 
which  encloses  the  target.  This  device  appears  adequate  for  laboratory 
experiments,  provided  the  expanding  target  effect  due  to  decreasing 
range  is  not  required.  This  sensor  has  been  used  to  design  and  evaluate 
the  performance  of  the  basic  servo  system,  since  basic  error  signals 
similar  to  the  CCT  tracker  output  are  available.  The  existing  gimbal 
system  is  driven  by  highly  geared  DC  motors  and  thus  the  typical  low 
frequency  problem  of  backlash,  stiction,  etc.  limits  performance.  A 
careful  quantitative  evaluation  of  the  tracking  accuracy  and  its  effect  in 
a  typical  homing  engagement  has  not  yet  been  performed;  however, 
some  measurements  have  been  made  of  tracking  performance  for  a 
circularly  moving  target  in  a  plane  perpendicular  to  the  sight  line.  The 
tracker  is  shown  in  Figure  20,  and  the  electronics  package  in  Figure  21. 
An  overall  view  of  the  open  loop  tracking  experiment  is  shown  in 
Figure  22. 
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In  the  design  of  any  servo  gimbal  system,  minimum  inertia  load 
is  desired.  Thus,  only  the  vidicon  and  its  deflection  coils  were  mounted 
initially.  However,  considerable  difficulty  was  encountered  in  trans¬ 
mitting  the  signal  to  the  remotely  located  preamplifier.  This  problem 
was  solved  by  mounting  a  preamplifier  stage  with  the  gimballed  vidicon 
assembly.  This  has  yielded  such  highly  satisfactory  results  that  future 
assemblies  will  probably  retain  this  feature. 

In  designing  the  servo  loop,  maximum  flexibility  was  gained  by 
using  commercial  DC  amplifiers  of  high  quality  for  signal  processing 
and  for  active  compensation  networks.  This  approach  was  selected 
because  the  dynamic  characteristics  of  the  CCT  tracker  are  not  yet 
known.  Loop  gains  and  compensation  network  characteristics  are 
easily  adjusted  in  the  experimental  equipment.  Furthermore,  if 
nonlinear  gain  functions  and  complex  compensation  become  desirable, 
these  are  easily  mechanized  with  the  operational  amplifiers. 


A  block  diagram  of  one  channel  of  the  servo  tracker  is  shown  in 
Figure  23.  The  sensor  transfer  function  (Walleye)  was  determined  to 
have  linear  DC  characteristics,  and  the  dynamic  response  could  be 
approximated  by  the  first  order  lag  shown  in  the  sensor  block.  Cancel¬ 
lation  compensation  was  employed,  and  the  lead  network  pole  was  placed 
as  far  away  as  possible  without  excessive  attenuation.  The  motor  and 
gear  train  transfer  function  were  derived  using  the  usual  approximations 
of  linear  torque -cur rent  characteristics  and  no  gear  backlash  or  stiction. 
The  numerical  transfer  function  was  obtained  from  the  formula: 


e_ 

V 


Kt 


4  ^(^  +  01 


where : 


Kt  =  Torque  constant 
Kn  =  Back  emf  constant 
Ra  =  Winding  resistance 
Jmz  Load  inertia 
N  =  Gear  ratio 


The  functional  schematic  diagram  of  the  servo  loop  is  shown  in 
Figure  24,  and  the  power  amplifier  schematic  in  Figure  25.  The 
current  output  versus  voltage  input  for  the  latter  is  shown  in  Figure  26. 
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Except  for  the  previously  mentioned  difficulty  with  slow  rates, 
performance  of  the  laboratory  system  seems  to  be  adequate.  Improve¬ 
ment  of  the  tracking  performance  in  the  marginal  region  by  pulse  width 
modulation  of  the  motor  driving  signal  has  been  investigated  but 
implementation  is  yet  incomplete. 

V.  CONCLUSIONS 

The  results  of  the  preliminary  analytical  and  experimental  work 
seem  to  indicate  that  the  principal  problem  areas  are  directly  related 
to  the  characteristics  of  the  TV  raster-scan  seeker.  Due  to  the  peculi¬ 
arities  of  this  device,  the  simple  concept  of  a  body-fixed,  fixed  field 
optics  tracker  does  not  seem  feasible.  At  the  minimum,  a  readily 
designed  gimballed  servo  mount  will  be  required. 

A  decisive  conclusion  cannot  yet  be  reached  on  the  terminal  guid¬ 
ance  problem.  It  seems  at  this  time  that  the  fielchof-view  problem 
cannot  be  resolved  within  the  basic  seeker  itself,  since  this  would 
require  an  increase  in  the  resolution  of  the  system. 

Future  planned  work  includes  careful  investigation  of  the  terminal 
guidance  problem  employing  the  dynamics  of  a  typical  missile  design, 
such  as  the  MAULER.  The  results  of  the  investigation  should  yield 
quantitative  data  on  miss  distance  as  a  function  of  coast  range.  If 
the  minimum  coast  range  required  by  the  seeker  does  not  result  in 
acceptable  terminal  performance,  and  state -of-art  refinements  of 
missile  design  (if  any  are  possible)  do  not  yield  sufficient  improvement 
then  the  only  alternative  would  seem  to  be  programmed  field-of-view 
optics. 

Studies  of  a  higher  performance,  compact  two-axis  servo  tracker 
have  begun.  Under  consideration  are  the  following  two  types: 

1.  A  device  employing  direct  drive  electric  torque  motors,  driven 
by  compact,  high  DC  amplifiers. 

2.  A  device  employing  high  performance  cold  gas  pneumatic 
actuators. 
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Appendix  A 


TEST  VEHICLE  CONSIDERATIONS 


A  survey  has  been  made  of  the  many  missiles  an-1  r’rones,  either 
in  existence  or  in  development,  to  determine  the  feasibility  of  using 
one  or  more  of  these  as  test  vehicles.  The  Ryan  Q-2  Firebee  drone 
seems  to  be  a  suitable  vehicle  for  open  and  perhaps  closed  loop  early 
flight  testing.  With  this  possibility  in  mind,  the  required  aerodynamics 
and  control  information  has  been  obtained  from  Ryan.  The  equations 
of  motion  for  this  type  vehicle  have  been  derived  and  are  shown,  along 
with  pertinent  numerical  data  and  graphs,  in  Appendix  B. 

Since  the  eventual  weight  and  volume  of  a  flight  configuration  of 
the  complete  CCT  seeker  cannot  be  accurately  estimated  at  this  time, 
the  selection  of  a  test  vehicle  for  testing  in  a  field  environment  becomes 
rather  difficult.  Therefore,  an  effort  has  been  initiated  to  obtain  all 
pertinent  information  on  several  missiles  in  a  range  of  sizes.  As 
these  data  are  obtained,  system  equations  and  simulator  programs 
will  be  written  for  each  missile  as  was  done  for  the  Q-2,  and  contact 
maintained  with  the  various  contractors  so  that  pertinent  changes  can 
be  incorporated  in  the  programs.  At  the  present  time,  effort  is  being 
concentrated  on  the  MAULER,  which  seems  to  offer  rather  good 
possibilities. 


Appendix.  B 


LINEARIZED  EQUATIONS 
OF  MOTION  APPLICABLE  TO  XQ2-C  DRONE 


Basic  Equations: 

1.  Force  ='mass  X  acceleration 

2.  Angular  acceleration  =  applied  moment/moment  of  inertia 

3.  Lift  (L)  on  an  airfoil:  L  -  */2  p  V2  S 
Drag  (D)  on  an  airfoil:  D  =  2/2  p  V2  S  Cq 

Where : 

CL,  Cjj  =  dimensionless  coefficients 
p  =  air  density 
V  =  air  speed 
S  =  airfoil  planform  area 


Assumptions: 

1.  Airplane  structure  is  rigid. 

2.  X,  Y,  Z  coordinates  rotate  with  airplane. 


3.  Disturbances  are  small  deviations  from  steady  state  flight. 

4.  Longitudinal  and  lateral  equations  are  decoupled. 

(See  Table  B-I  for  symbols  and  definitions) 


Using  Newton's  Laws: 

0  sin  (p 
0  cos  tfi  . 

•k  *“  —  ^xz*" 

M  =  IyP 

N  =  I  r  -  I  P 
1zr  xz 


m(U)  =  X  -  mg  sin  0 
m(V  +  Vr  =  Y  +  mg  cos 
m(<l>  -  Vq  =  Z  +  mg  cos 


I 


II 
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Equation  I  describes  translational  motions,  whereas  Equation  II  denotes 
angular  motions. 


Table  B-I 

DEFINITIONS  AND  SYMBOLS 
Symbol  Definition 


X,  Y,  Z 

Forces  along  X,  Y,  Z  axis 

L,  M,  N 

Moments  about  X,  Y,  Z  axis 

U,  V,  W 

Flight  path  velocities  along  X,  Y, 

Z  axis 

p,  q,  s 

Components  of  angular  velocity 
about  X,  Y,  Z  axis 

e 

Pitch  attitude  angle 

Roll  angle 

Yaw  angle 

V 

Flight  path  angle 

a 

Angle  of  attack 

*a 

Aileron  deflection 

Elevator  deflection 

6r 

Rudder  deflection 

T 

Aircraft  time  constant 

Tv 

Change  in  thrust/unit  speed  change 

V 

True  air  speed 

b 

Wing  span 

C 

Mean  aerodynamic  chord 

^xx’  *yy'  *zz 

Momenta  of  inertia  about  X,  Y,  Z 
axis 

Kx>  Ky>  Kz 

Radius  of  gyration  about  X,  Y,  Z 
axis 

Cij 

Aerodynamic  coefficients  where  first 

subscript  (i)  denotes  type  of  forcing 
function  and  the  second  (j)  denotes 
the  angle  or  angular  rate  associated 
with  the  coefficient 


Table  B-I.  (Concluded) 


Symbol 
i  =  m 
i  =  n 
i  =  L 
i  =  l 

i  =  y 


Definition 

Pitching  moment  coefficients 
Rolling  moment  coefficients 
Lift  coefficients 
Yawing  moment  coefficients 
Side  force  coefficients 


If  one  assumes  small  disturbances  from  steady  state  flight,,  the  longi¬ 
tudinal  equations  of  motion  of  an  airplane  can  be  written  in  terms  of 
the  stability  derivatives: 


Lift  Equation  ( sF  =  ma) 


mZ  =  V2  p  V2  S  Cha  a  +  2  (%p  V2  SCLU). 

Using  the  approximation 

mZ  =  mVy  ,  mVy  =  %  p  V2  S  +  2  (%  p  V2  S  CL  U)  . 

Above  equation  is  more  easily  applied  if  rewritten  explicitly  as  a 
function  of  the  aerodynamic  coefficients. 

Using  the  identity  9  =  y+  &  or  y  =  6  -  a 

m V(fl  -  a)  =  %  p  V2  S  Chaot  +  2  (%  p  V2  S  CLU)  . 


Transposing  the  left  side  of  equation  to  the  right: 

(Vi  P  V2  S  Clo,®)  +  2  (Vi  p  V2  S  CL  U)  +  mVd  -m  V0  =  0  . 
Multiplying  through  by  2  and  dividing  by  p  S  V 


vcl*«  +  2V  clu  =  e 


Dividing  through  by  -V; 


_CL0'Q' 


-  2Clu  - 


2m 


■  a  + 


2m 


pVS  pVS 


6  =  0 


The  term  —77-7  has  the  unit  of  seconds  and  can  be  defined  as  the  air 
pVS 

aircraft  time  constant  (t). 

Thus  the  lift  equation  becomes: 
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-cLaa-  ZCL  U  -  Zi'ot  +  2t6  =0.  (A) 

Moment  Equation  (2M  =10) 

V2P  V2  S  C  (Cmi  CmSe4a)  =  mkx2  6 

Multiplying  through  by  2  and  dividing  by  p  V2  S  C  the  following  is  obtained: 

"m.0  0  +  ^m6e 

where  is  again  the  time  constant. 

Thus  the  moment  equation  can  be  written: 


cmi  cm5e  “  ~  (^)!>  =  0  ■ 


mot 


2m 

pVS 


m 


Z  •  or  +  C  ot  —  C  vfl  +  C  p  6e  —  4t  6 

ma  mo  m0  m6e 


I7Xt)  »•  ■  <*> 


The  longitudinal  flight  path  equation  can  be  written: 

V"[mU  +  Tv  u]  =  2(VaP  V2  S)<CD  U  +  CLV) 

-  Y2P  V2  S  (Clq.-C-q^  -CDq)  • 

2  m 

Multiplying  by  and  substituting  for  —  — g  : 

-2tU+^  Tv  U  -  2Cd  V  -  2ClU  +  C^Q2  -  CD„ 


~  cLej0  =  °- 


(C) 


The  lateral  equations  can  be  obtained  in  the  same  manner  and  are 
given  in  our  standard  form: 

-  2t  p  +  Cyp  p  +  Cj_,<ji>  -  2t  4<  =  0  (D) 

C/P  o  -  4^XY)!  *  +  Ci^  +  Ca+  +  C»a  =  0  <E> 

C„p  P  +  ^  <J>  -  4t(£)(*s)‘  Cn6r  Sr  =  0  (F) 

Equations  A  through  F  can  then  be  programmed  on  the  analog  computer 
and  are,  when  properly  scaled,  valid  for  the  simulation  of  the  Q2C 
aircraft  within  the  limitations  imposed  by  initial  assumptions.  However, 
before  the  above  equations  can  be  conveniently  programmed,  other 
assumptions  need  be  made,  particularly  with  respect  to  the  stability 
derivatives:  Many  of  the  derivatives  are  functions  of  several  variables, 
such  as  Mach  number,  center-of-gravity  station,  angle  of  attack,  etc. 
However,  most  of  the  derivatives  become  independent,  or  at  least 


23 


linear  functions  of  speed  if  one  restricts  himself  to  speeds  below 
Mach  0.9.  One  can  further  simplify  the  problem  by  assuming  airspeed 
to.be  a  constant,  since  most  derivatives  are  then  reduced  to  constants. 
If  more  than  one  speed  is  of  interest,  several  independent  studies  can 
be  made  for  incremental  values  of  airspeed.  Assumption  must  also 
be  made  for  altitude  and  payload,  since  these  influence  the  aircraft 
time  constant  (due  to  p)  and  center  of  gravity  station. 

Pertinent  data  for  the  Q2-C  are  shown  for  three  velocities 
(Mach  0.6,  0.8,  and  0.9)  in  Table  B-II.  English  units  are  used 
throughout. 


Table  B-II 

Q-2  NUMERICAL  PARAMETERS 


Altitude 

V  =  Mach  0.  6 
Sea  level 

V  =  Mach  0.  8 
Sea  level 

V  =  Mach  0.  9 
5,  000  ft 

Gross  weight 

2, 060  lb 

2,  060  lb 

2,  060  lb 

Mass,  slugs 

64 

64 

64 

P 

0.00237 

0. 00237 

0.  0020 

b 

12.  9  ft 

12.9  ft 

12.  9  ft 

S 

39.  9  ft2 

39.  9  ft2 

39.  9  ft2 

V 

690  ft/ sec 

920  ft/sec 

1 ,  050  ft/ sec 

r 

1  sec 

I  sec 

0.  85  sec 

C  station 
g 

88. 5  in 

88.  5  in 

88. 5  in 

O' 

0 

0.  5  deg 

0.75  deg 

0.  7  deg 

c 

2.783  ft 

2.  783  ft 

2.  783  ft 

cT  t  . 

0.  067 

0.074 

0.  073 

La?  trim 

^Da? 

See  Cq  versus 

a?  curves 

T 

A  V 

Negligible  for  jet  aircraft 

CDo 

0.  002 

0.  002 

0 

cLa 

0. 08/deg 

0.  088/deg 

0. 085/deg 

Cma? 

-0. 034/deg 

-0. 032/deg 

-0. 032/deg 

r 

-4. 6/rad 

-7/rad 

-4. 2/ rad 

cm6 

-17/rad 

-20 . 1/rad 

-16. 2/ rad 
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Table  B-II  (Goncluded) 


Altitude 

V  =  Mach  0.  6 

Sea  level 

V  =  Mach  0.  8 

Sea  level 

V  =  Mach  0.  9 

5,  000  ft 

°m6e 

-02/ deg 

-0. 0181/ deg 

-0. 0165/ deg 

GyP 

-0.  0225/deg 

-0. 0245/deg 

-0. 0265/deg 

Cl(p 

-0. 37/rad 

-0. 037/rad 

-0. 037/rad 

CZ6a 

0.  0005/deg 

0. 00057/deg 

0. 00055/deg 

C  np 

0.  0035/deg 

0. 0038/deg 

0. 0039/deg 

Cnijj 

-0. 63/ rad 

-0. 63/ rad 

-0.73/rad 

Cn6r 

-0.  00077/deg 

-0. 00075/deg 

-0. 00075/deg 

Kx 

0.987  ft 

0.  987  ft 

Ky 

3.89  ft 

3.  89  ft 

Kz 

3.87  ft 

3.  87  ft 

62  slug/ft2 

62  slug/ft2 

h 

1 . 036  slug/ft2 

1.  036  slug/ft2 

1 . 025  slug/ft2 

1.025  slug/ft2 

^xz 

46.  5  slug/ft2 

46.  5  slug/ft2 

The  differential  equation  form  is  most  suitable  for  analog  pro¬ 
gramming;  however,  for  checking  and  analysis,  the  transfer  function 
form  is  desirable.  Derivation  of  6/6e( S),  <£/6a(S)  i|j/6r(S)  follows: 
Rewriting  the  longitudinal  equations  in  operation  form,  using  S  =  d/dt: 

(-C;lq/-2t  S)  of  -  2Cj_U  +  2t  S  9  -  0  (A1) 

(^ma  S  +  CmQ,)  ot  +  (-C S  —  \  S2)  9  =  -Cm6e  6e  (B1) 

where 

>  ■  nmf 

[(t)  -  2t  S)  -2Cq  -2Gl]  U  +  (C^  -  C-q^)  or  -  -  0  (C1) 

where 

*  m  v 

The  determinant  solution  for  6/ 6e  is  given  by: 
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\ 

] 

I 


_(S)  = 
6e 


(-  CLa  -  2t  S)  -  2Cl  0 

(^mi  ®  +  ^mo?)  ®  _<-'m6e 

(cLa  ~  cDcy)  (g  -  2t  S  -  2Cp  -  2CL)  0 


w'he  r  e 


V  =  4t  I 


Cnp  P  +  Cn<^S  <t>  -  (V  S2  -  cn^i  S)  4»  =  -  Cn6r  6r 


where 


X."  =  4r 


b\/K'2 


The  determinant  solution  for  (f)  (S)  is  then  given  by: 

(-  2r  S  +  Cyp)  0  -  2t  S 

cLp  ~Cl6a.  6a  C24<  S 

Cnp  -Cn6r  Sr  -(V  S2  -  Cmf  S) 


W  = 


(-  Zt  S  +  cyp)  ct  -  2r  S 

-(V  s2  -  c l(p  S)  c 4  s 


cLp 


Cnp 

Similarly,  t)j  (S)  becomes: 

(-  2t  S  +  Cyp) 


+  (S)  = 


<-'n<£  S 


-(V  S2  -  Cn^S) 


c/p 

Cnp 


-{V  S2  -  C l(p  S) 
Cn0  S 


~  CZ6a  6a 
—  Cn6r  Sr 


-  (2t  S  +  Cyp)  - 


C  l 


2t  S 


C/^  S 


Cjp  -(V  S2  -  Cl(p  S) 

Cnp  Cn^  S  -<X"  S2  -  Cm [  S) 


(-  Cl-  -  2t  S)  -  2Cl  2t  S 

(CmiS+CmCI)  0  (~  Cm0  s  -  xs2)| 

(cLa  -  cDa)  (t|  -  2r  S  -  2Cq  -  2CjJ  ~  CL 
The  lateral  equations  can  similarly  be  written: 

(-2 t  S  +  Cyp)  p  +  Cj0  -  2t  S  i)j  =  0 

czp  P  -  (X-  S2  -  C^  S)  4>  +  C^  S  *  =  -  Cz6a  6a 
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The  equations  derived  herein  are  sufficient  for  aerodynamics 
simulation  of  unsymmetric  (bank  to  turn)  vehicles.  In  the  case  of  the 
Q2-C  drone,  all  of  the  required  numerical  values  are  available.  For 
a  study  of  the  complete  loop,  including  the  seeker,  the  Q2-C  autopilot 
transfer  function  is  required.  These  additional  data  are  given  as 
follows: 


Elevator  inner  loop: 


(l  +  il)(1  +To+J§w) 


Aileron  inner  loop: 


(100W 


i 


j)  ‘jsS-re;  Aq  papuajqns  sjSuv 
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Q  <50§j'e;  Aq  papuaijqns  ©{Sup  jo  aStreqo  jo  0^-b>j 
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Figure  4.  PLOT  OF  df  VE 
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igure  8.  BLOCK  DIAGRAM  OF  SIMPLE  HOMING  LOOP 


Figure  9.  COMPUTE 
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Figure  20.  PHOTOGRAPH  SHOWING  TV  SERVO  TRACKER. 
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Figure  21.  PHOTOGRAPH  SHOWING  TRACKER  ELECTRONICS. 


Figure  22.  PHOTOGRAPH  SHOWING  TRACKING  EXPERIMENT. 


-25  volts 
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Figure  25.  MOTOR  DRIVER 
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Figure  B-3.  DRAG  VERSUS  ANGLE  OF  ATTACK. 
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